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Abstract 

Using the Eliashberg strong coupling theory with vertex correction, we calculate 
maps of transition temperatures ( Tc ) of electron-phonon superconductors in full 
parameter space. The maximums of transition temperatures for superconductors are 
predicted based on the maps and the criterion of instability of superconductivity. 
The strong vertex correction and high transition temperature are tightly correlated 
in superconductors. We predict that the maximum of Tc of new iron-based super- 
conductors will be close to 90 (K). 
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1 Introductions 



The researches to find home-temperature superconductors are still highlight 
in the field of material science [1]. The theoretical predictions of Tc are the 
main efforts of theoretical scientists working in this field. The BCS theory [2] 
and Eliashberg strong coupling theory [3-7] are well known acceptable the- 
ories to explain the superconductivity of electron-phonon superconductors. 
There exists the well known upper limit of Tc defined by the McMillan Tc for- 
mula [7]. The lattice instability (softening phonon) and the total strength of 
electron-phonon coupling also set the bounds of Tc of electron-phonon super- 
conductors [8]. The dielectric properties have been expected to have significant 
influence on the possible maximum of Tc for a variety of superconductors, but 
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the problem is still not answered because the complex relation between ef- 
fective electron-electron interaction and dielectric function [9]. The magnetic 
fluctuation and the vertex correction (non-adiabatic effects) [10] can suppress 
the Tc and lead to the instabihty of superconductivity of electron-phonon 
superconductors. . 

The new novel Iron-based high-temperature superconductor [11-15] has been 
discovered and the new record of transition temperature Tc has reached about 
57.4 (K) by doping rare earth elements into CaFeAsF [15]. Especially the sig- 
nificant normal isotope effects of iron with a ~ 0.4 for Bai_a;K^Fe2As2 indi- 
cates that electron-phonon interaction plays an essential role to explain their 
superconductivity [16]. The Eliashberg theory suitable for both weak-coupling 
and strong-coupling should be studied in more detailed. It is very desirable to 
study the effects of vertex correction beyond Migdal theorem [10], especially 
for superconductors have components of light atoms. The Eliashberg theory 
has been successfully used to calculate Tc of many types of superconduc- 
tors [17-19]. However, only special regions in parameter space of A — f2p — /x* 
have been explored, where A is the parameter of electron-phonon interaction, 
Qp the frequency of phonon and //* the Coulomb pseudo-potential. In this pa- 
per, we present the information ( Tc map ) in full parameter space A — fip — 
The Tc maps obtained in this paper are very helpful to analyze the relation 
between Tc and these superconducting parameters and guide to find the su- 
perconductors with higher Tg. 



2 Theory 

In the paper, we generalize the equation of energy gap in reference [10] by 
including the Coulomb interaction. The standard energy-gap equation with 
the form of imaginary Matsubara's formulation [5] and the equation after 
having considered vertex corrections are used in our calculations. For isotropic 
electron-phonon interaction. 



the calculations of vertex corrections are greatly simplified. The electron- 
phonon interactions are included in the vertex corrections only by the func- 
tions of electron-phonon interaction A(n) defined below. The self-energy and 
Green's function of electron in the Nambu scheme are expressed as 
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where (f){k,iujn) is the function of energy gap, Z{k,iuJn) the renormahzation 
function, and SkiiUn) = eKiUn) + xik^iuj^). In above equations fi and are 
two Pauh matrixes and tq 2x2 unity matrix. 

If the density of state of phonon F{v) = V/{2ttY J d'^pB{p, u) and its spectral 
function B{p, v) are known, the effective electron-phonon interaction a^(z/) is 
defined by 



c^»F{u) = ^[^[ ^B{p-p',u)\j;^^,\\ (3) 
W Js vp JS' vp' 

where W = {2ttY Jg d'^'^p/vp, Fermi velocity and </p_p' the matrix element 
of electron-phonon interaction which is constant J in this work. The first 
order self-energies of electrons contributed from electron-phonon interaction 
and electron-electron interaction are standard and written as 

rEp 

T.\iuj„) = kBT^\{n-n') j dek'hG^ek'.iuJn')^ (4) 
T,''{iujn)= -ksTj^l^* / dek'f'iG°'^{ek',iuJn')% 



where G°'^{ek', ioOn') is off-diagonal part of Green's function matrix G{eki, iuJn')- 
For convenience, the kinetic energy e^/ has to shift Sk' — Ep and a symmet- 
ric band from —Ep to Ep is chosen. Additionally we assume /c'-dependent 
G{k', iun') realized by G{ek', iuJn')- The electron-phonon interaction is included 
in X{n) = 2 /q°° (iz/a^F(z/)z//(z/^ + cj^). The Coulomb pseudo-potential is de- 
fined as /i* = Ho/ {1 + fio\n{Ep / uo)) , where fiQ = N{Ep)U, N{Ep) the density 
of state at Fermi energy Ep,l] the Coulomb interaction between electrons and 
Uq characteristic energy of typical phonon correlated to superconductivity. 

The lowest order vertex correction comes from the second order self-energy 
that is written as 



T,'^{k,iuJn) = {kBTf\J\^T.n'n"k'k"D{k - k' ,iVn~n')D{k - k",Wn-n") (5) 

X %G{k', iujn')hG{k"' , iuJri"')hG{k" , iuJn")h 

where D{k* ,ivn*) is phonon Green's function with the definitions of indexes 
n*=n — n' or — n", n'" = n' + n" — ra, k*=k — k' ot k — k" and k'" = 
k' + k" — k. The summation J2k* of wave vector k* {k* = k' or k") is re- 
placed by V/{27i)^ J d^k* = V/vp{27i)^ J d'^k* J dsk*. We introduce an integral 
{l/2Ep) J^p^ dsk'" in Eq.Q. So the second order self-energy after having av- 
eraged on Fermi surface can be simplified as 
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= ^-^^^n'n"X{n - n')X{n - n") (6) 
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now Ski" is an independent variable. Three integrals of kinetic energy can 
independently be done by 

^ An 

dek*G{ek*,iuJn) = -7r(is„ro + -p^ri)a„ (7) 

-E* \UJr,.\ 



when temperature is very close to T^, (l){k,iujn) — > 0, so we ignore the term 
iUn) in the denominator of Green's function Eq.([2]). Other parameters are 
defined as s„ = Un/\uJn\, o-n = (2/7r) arctan(£'*/|co'„|) and A„ = (piiojn) I Z{iujn) 
the energy-gap parameter. In the general situation, the second-order self- 
energy cannot simply be written as the functions of \{n). In terms of above 
approximation, we transform the formulas of second-order self energy to a 
new formation which is suitable for numerical calculations. Although effects 
of electron-electron correlation are not treated accurately, the above approxi- 
mation is still reasonable because the main contributions to superconductivity 
come from \{n — n')\{n — n") in Eq.(l5]) and Eq.(l6]). 

The total self-energy we consider is 'L{iuJn) = T}{iujn) + S'^(ic<j„) + T?{iun). 
By combining with Dyson Equation and considering that, when temperature 
is very close to Tc, A„ 0, the terms proportional to A^ are ignored and the 
energy-gap equation is generalized to 

{Knn' - pSnn')-r^ = 0) (§) 

, \UJr,l 



where the kernel matrix 



Knn' = [A(n - n')B{nn') - fi* + G{nn')]an' - Snn'Hn', (9) 
+00 1^ ^1 

Hn'= Y] [6n'n"—r-7^ + Hn' -n")A{n'n")sn'Sn"an"] 

n"=—oo ^ 

with parameters A{nn') = 1 — VA(nn'), B{nn') = 1 — VB(nn') and a„ = 
(2/7r) arctan(ii^^/Z„|a;„|). The 3-dimensional parameters of vertex correction 
have the form 
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Vsinn') 



VA{nn') 



'^ksT /2Ef) ^ \{n — n")Sn'+n"~nSn"(ln'+n"-nO". 
n" 

"^ksT / Ep) ^ X{n — n") Sn' +n" -nSn" O^n' +n" -nO'n' 



•n' 



It 



(10) 



C{nn') 



^kBT/2EF) X{n - n")X{n' - n")sn'-n"+nS, 




n" 



If we ignore the vertex corrections, three parameters V^(rm'), VB{nn') and 
C{nn') are all equal to zero and the kernel Eq.Q of energy-gap equation 
reduces to the general form without vertex correction [5]. At the same time, the 
differences of the equations of energy gap in 3-dimensional and 2-dimensional 
space are only included in the function a^F^v) if vertex corrections are ignored. 
The introduction of pair-breaking parameter p creates an eigenvalue problem. 
The physical gap-equation is corresponding to p=0. In the calculation of a„, we 
choose Zn ~1 the value of normal state. The T^ is defined as the temperature 
when the maximum of eigenvalues E'""'^' of kernel matrix Knn' crosses zero and 
changes its sign. We use about N=200 Matsubara's energies to solve above 
equation. The method to find Tc is similar to the bisection method used to 
find the roots of a nonlinear equation [20]. The nonlinear equation to find 
Tc is E™"^ (T)=0. We can reach the accurate of 0.0001 (K) after only 20-30 
iterations in temperature interval from (K) to 600 (K). 

In calculations of a^F(z/), we assume that a^{v) is approximately a constant 
around the peak of phonon mode and the density of state of phonon is ex- 
pressed as 



where f2p is the energy of phonon mode, 002 the half-width of peak of phonon 
mode and uj^ = 2uj2- The parameter of electron-phonon interaction is defined 
as A = A(0) = 2 dva'^F{v)/v. An important relation 



will be used, where M is the effective mass for a certain phonon mode, espe- 
cially the constant r] , the so-called McMillan-Hopfiled parameter, is closely 
related to Tc [5]. The superconductor parameter rj = N{Ep){J'^), character- 
izes the chemical environment of atoms and almost keeps as a constant against 
the simple structural changes or isotope substitutions. 




(11) 



M{u^)\ = 1] = const 



(12) 
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Fig. 1. The 3-dimensional Tc map and three slice planes with A=2.0, /u*=0.1 and 
i7p=80 meV. The numbers within the rectangles nearby the contour lines are Tc 
values of the corresponding contour lines. 

3 Tc Map without vertex corrections 

A Tc map in full parameter-space X — Qp — fi* is very helpful to understand 
the superconductivity with electron-phonon mechanism. At first, we choose 
three-dimensional 30x30x30 mesh-grid with 0.4< A <4.0, 0< fi* <0.4 and 
5.0 meV< Qp <80 meV. We calculate Tc on every mesh point using Eq.([S]). 
We plot three slice planes with A=2.0, yU*=0.1 and flp=SO meV respectively 
shown in FiglH Comparing with others two parameters, the choice of /i* has 
smaller effects on Tc. The parameter A=2.0 is the upper limit of electron- 
phonon interaction [7]. When the parameter A is larger than 2.0, the crystal 
lattice will be unstable. In our calculations, the half-height width uj2 has the 
same value (4 meV ) for all points of the Tc map. If the peak of phonon mode is 
not too broad, the Tc map has no significant change. Conversely, the Einstein 
mode with 002 ~0 has higher Tc. 

In Figl^l we plot a Tc map in two-dimensional X — Qp plane in details using 
50x100 mesh-grid with parameters A from 0.4 to 4.0, Qp from 5 meV to 160 
meV and the Coulomb pseudo-potential n*=0.1. The two-dimensional map 
is the extension of the slice of the corresponding 3-dimensional map with 
/i*=0.1 shown in Fig{Tl The important role of the relation M(z/^)A=const 
of Eq. flT^ can be shown clearly on the Tc map. If the effective mass M of 
a certain phonon mode is constant, so is the parameter {y'^)X. We assume 
the half-width of peak of phonon mode uj2 -C VLp, so {y'^)X ~ f2pA=const. 
The two curves of f2pA=320 and 3600 (meV)^ are plotted in the same map. 



In the regime of strong coupling, ksTc is approximately equal to aJ {u'^)X 
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Fig. 2. The Tc map in X — Qp plane with //*=0.10 obtained from imaginary-energy 



~ ayfipA, where a=0.115. Thus (z/^)A=const is equivalent to Tc=const. In 

fact, (z/^)A=const was used to define the possible maximum of Tc [7]. The 
Tc increases with A and saturates at A = 2 by using the McMillan formula. 
Similar to McMillan's idea, our results in the Figj2]just show that the two 
curves coincide approximately with the contour lines of Tc in strong-coupling 
regime. The approximate saturations happen when A is very close to 2. 

The term (z/^) A characterizes the possible maximum of Tc for fixed (z^^)A value 
if we consider possible bounds on A. We can see from the FigJ21that the curve 
will have more chances to meet contour lines with higher Tc for larger 
Thus, the maximums of Tc above two curves with f2pA=320, 3600 (meV)^ 
are close to 20K and 80K respectively. When superconductor under small 
structural modifications or isotope substitution, the corresponding parameters 
A and VLp move along these curves because the chemical environments keep 
almost unchanged. In the region of high-energy phonon, the smaller change of 
A can induce larger change of Tc because the curve npA=3600 (meV)^ spans 
more contour-lines. 

To study effects of Coulomb pseudo-potential /i*, we calculate two others Tc 
maps in /i* — A plane with mesh-grid 50x50 for two phonon energies f2p=10 
meV and 40 meV respectively (Fig 13]). For f2p=10 meV, the maximum of Tc 
is about 24 K, however 100 K for i7p=40 meV. We can find that when /x* > 
0.2, II* have smaller infiuences on Tc. 
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Fig. 4. The effects of vertex correlations on Tc map (//*=0.10) with two parameters 
c=0.12(eV)/Ei;'=0.03 and 0.06 compared with the results without vertex corrections 
c=0.0 (Ep = oo). Two groups of contour lines with Tc=70K and 140K are ploted 
for three values of c=0.0, 0.03 and 0.06 respectively. The inserted figure plots the 
curve of Tc by c at a fixed point P shown on the X — 0,p plane. 

4 Tc Map with vertex corrections 



It is very important to study tlie vertex correction in general strong coupling 
theory [10]. We control the vertex correction by changing the Fermi-energy 
Ep in the parameter a„ = (2/7r) arctan(£^i?/Z„|a;„|). We choose two values 
1.0 eV and 2.0 eV of Ep. If Ep is too large 2 eV) the vertex corrections 
are invisible in our calculations. Two groups of contour lines with Tc=70 K 
and 140 K are plotted, and for each of groups, there are three contour lines 
for c=0.0, 0.03 and 0.06 respectively, where c=0.12(eV)/Eir is the control 
parameter of vertex correction, especially c=0 or E^? = oo is corresponding 
to the absence of vertex corrections. We can see that for two Tc, the vertex 
corrections are so small that the Tc map has no significant changes. For the 
same Tc, the contour lines only move to regions with larger A slightly. In the 
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Fig. 5. The positive (A >3.2) and negative vertex corrections (A <3.2) at high 
phonon energy O,p=90 meV. 

inserted figure of FigJH at a point P witli fixed A=1.50 and f2p=65meV, the 
vertex corrections decrease Tc to smaller value. However, the negative vertex 
corrections only happen when Qp <80 meV and A <4. Beyond the ranges of 
parameters, such as Qp=90 meV and A >3.2, the positive vertex corrections 
(Tc increases with uo/Ef) are realized (Figl5]). It's not surprised that, if A <3.2, 
only the general negative vertex-corrections are found although the phonon 
energy Qp is larger that 80 meV (Figj5]). The key point is that, although the 
effects of vertex corrections are rather complicated, they do not change the T^, 
maps significantly if Ep isn't too small. 

We expect that the vertex corrections are strong for superconductors with 
high Tc. From the FiglH we can see that, when f2p is smaller than 40 meV, 
the vertex effect is very small although the parameter A of electron-phonon 
interaction is probably large. The effect of vertex correction is determined by 
(z/^)A and not solely by the parameter A of electron-phonon interaction in 
the region of low phonon-energy. The vertex correction is closely correlated 
to Tc because of Tc oc y^(i^^)A in strong coupling regime with A > 2.0. We 
will prove that the vertex correction is also characterized by (i^^)A in the 
region of low phonon energy. It easily sees that the parameter VAinn') = 

Pv J2m X{l^m)Sn'-mSn-man'-man-m wherC Py = '7l'^kBTj2EF and 171 = 12- u" . 

If the Einstein spectrum a^F(z/) = (A/2)f2p5(z/ — Qp) is considered, we can 

get 

X{iy„^) = \nl/iul + nl). (13) 

If Qp 0, then X^Urn) {^/^m)^^py the parameter of vertex correction 
V{nn') PvC{n, n', Tc)Afi|„ where C{n, n', T^) = Y.m Sn'-mSn-man'-man-m./ 
is well defined and finite. So the parameter f2pA or (i^^)A characterizes the ver- 
tex correction. In real materials f2p(<C Bp) has an up limit Vtp,- When Vtp ap- 
proaching to high energy, Xiym) -^^D/l^^m+^i)); yip^n') PvD{n,n' ,Tc)X, 
where D{n,n',Tc) = J^mSn'-mSn-man'-man-rnQD/i^m + ^d) Well defined. 
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So we can see that the vertex correction is characterized by A at high energy. 
From the Fig|2l is almost independent on Qp when f2p > 40 meV and 
A < 2.0, and is determined by A as well. Thus, we have proven that the 
vertex correction (as well as Tc) is characterized by (z^^)A in the region of low 
phonon energy and by A in the region of high phonon energy. There is a char- 
acteristic energy close to 40 meV and the vertex corrections (as well as Tc) 
have different behavior when fip is larger or smaller than the characteristic 
energy. 



5 Discussion and Summary 

For superconductors with Tc close to home temperature, the parameter A of 
electron-phonon interaction should be larger than 2.0 or the phonon-energy Qp 
larger than 80 meV (Fig. [2]). Generally, we hope to find higher Tc by increasing 
the frequency of phonon by choosing compounds including light atoms such as 
Hydrogen atoms in molecule crystals [21]. The contour line Tc=300K in Fig 12] 
is beneath the line of A =2.0 when Qp is larger than 150 meV or 1210 cm~^ 
for Raman shift. We can see that only small changes of A can induce large 
changes of Tc if fip > 80 meV. In compounds with components of light atoms 
and heavy atoms, the frequencies of optical modes are determined by light 
atoms. The heavy atoms have the capability to stabilize the lattice vibrations 
so it is advantaged to form macroscopic-coherent superconducting states. 

It is very interesting that we can explain the spatial anti-correlation between 
energy-gap and phonon energy of cuprate superconductor Bi2212 [22] in terms 
of the electron-phonon mechanism by having solved the real-energy Eliash- 
berg Equation under the constraint from Eq. f[T2]) [23]. The large value of A 
(>1.0) can be obtained from the frozen-phonon method by considering non- 
local long-range Madelung interaction which is non-adiabatic effect when an 
atom displaces away from its equilibrium position [24,25]. Our results in this 
context also indicate that the vertex corrections which include part of non- 
adiabatic effect are closely correlated to high Tc. On the other hand, in terms 
of the linear temperature-dependent resistivity at high temperature above Tc, 
the strong electron-phonon coupling could be extracted by the analysis of 
strong coupling theory [26,27]. The upper limit of phonon energies of in-plane 
breathing modes of cuprate superconductors is close to f2p=80 meV [28]. If 
the anisotropy of Bi2212 crystal is ignored, from the Fig{Tl we can see that the 
maximum of Tc will be within 140K-160K. The highest Tc is about 130 K for 
cuprate superconductors close to the range. From the FigJU there are larger 
vertex corrections at higher Qp and larger A, however the Tc doesn't changes 
significantly because the band-width of Cu02 plane for cuprate superconduc- 
tors such as Bi2212 is not too small about 4 eV, and the Fermi-energy Ep is 
at least 1 eV [29,30]. 
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Finally, we discuss the new iron-based layered high-temperature supercon- 
ductors [11-15]. The carriers are mostly possible in FeAs layers. The multi- 
band electronic structure, especially the cooper-pairs spanning different Fermi- 
surface sheets probably play an important role to its superconductivity. The 
pairing mechanism of different Fermi-surface sheets is absent in present model, 
but the model is still equivalent to the mean of multi-band model. Very simi- 
lar to the cuprate superconductors, the linear-response calculation of electron- 
phonon interaction only provides very small parameter A of electron-phonon 
coupling that does not account the requirement of high T,, from 27(K) to 
57(K) [31]. However, just like for cuprate superconductors [25], the nonlocal 
long-range Madelung interaction generally enhances parameter A of electron- 
phonon interaction and account the requirements for high Tc. Experimental 
evidence for strong electron-phonon coupling comes from the measurement of 
normal-state resistivity at high temperature [32]. The upper limit of phonon- 
energy is about 40 meV [33]. From the Figl2]and Fig|31^b), we can predict that 
the maximum of transition temperature Tc for layered iron-based supercon- 
ductors may access to 90 (K). 

In summary, we calculate Tc maps in full parameter space using strong cou- 
pling theory including vertex corrections. We also discuss the possibility in- 
creasing Tc by using materials including light atoms with higher phonon fre- 
quency. The strong vertex corrections are inevitable in materials with high Tc. 
Based on the maps obtained in this work, we predict the maximum of Tc of 
new finding iron-based superconductors is about 90 (K). 
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